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Abstract. Quantum computation in the one-way model requires the prepara- 
tion of certain resource states known as cluster states. We describe how the 
construction of continuous-variable cluster states for optical quantum computing 
relate to the existence of certain families of matrices. The relevant matrices are 
known as weighing matrices, with a few additional constraints. We prove some 
results regarding the structure of these matrices, and their associated graphs. 



PACS numbers: 02.10.Ox, 03.67.Lx, 42.50.Ex 

1. Introduction 

In the standard model of quantum computing [1], the computation proceeds by 
coherently performing unitary dynamics on a simple initial state of a quantum system 
before being read out by a sequence of local measurements. By contrast, the one- 
way model of quantum computing [2] eliminates the need for coherent quantum 
evolution, substituting instead a highly entangled, but easily prepared initial state 
as the resource for the computation. The typical resource state used is the cluster 
state [3], whose preparation consists of, for example, coupling a system of qubits with 
a nearest-neighbor Ising interaction. The computation proceeds by a sequence of local 
measurements on the cluster state and the basis chosen at each step determines the 
computation. In this way, a cluster state is like a quantum breadboard on which the 
quantum circuit is inscribed via measurements. 

Although cluster states were originally defined using qubits, they generalize to d- 
level quantum systems [4] and continuous variables [5, 6]. Recently, Ref. [7] proposed 
a method for efficiently generating continuous- variable cluster states (CVCS). The 
method relies on pumping an optical cavity containing a special nonlinear medium at 
certain carefully chosen frequencies. The resulting CVCS is encoded in the quadratures 
of the photons that emerge from the cavity. 

The method of CVCS generation proposed in [7] required, on grounds of 
experimental and theoretical tractability, the existence of certain families of matrices. 
The required matrices are weighing matrices, defined as matrices W whose entries 
are in {0, ±1} and satisfy WW T = kl, where k is called the weight. The weighing 
matrices were obliged to have Hankel form, meaning that the skew-diagonals of the 
matrix were all constant, and any skew-diagonal with support on the main diagonal 
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must be identically 0. A matrix with vanishing main diagonal is said to be hollow. 
As we show later, when viewing the matrices as adjacency matrices for an associated 
graph, this is implied when the graph is connected and bipartite. 

In this article, we investigate the mathematical structure behind the construction 
in [7] and begin to classify the associated matrices. The discussion is organized as 
follows. In section 2 we give a brief glimpse of the physics which motivates the 
restrictions given in the previous paragraph. We begin the mathematical discussion 
in section 3, where we discuss anticirculant graphs, and section 4, where we discuss 
anticirculant weighing matrices, a special case of the Hankel weighing matrices. In 
section 5, we complete the rather trivial classification of the case of weight 2, while 
section 6 provides several examples about the case of weight 4. We discuss a number 
of open problems in section 7, and conclude in section 8. 

2. Quantum computing in the frequency comb 

In this section we outline the physics that naturally leads us to consider Hankel hollow 
weighing matrices. For a more detailed discussion of the physics, see [7, 8, 9, 10]. 

The resonant frequencies of an optical cavity arc defined classically as the modes 
which constructively interfere inside the cavity. Because they are evenly spaced, they 
form a so-called a frequency comb. By placing a nonlinear medium inside the cavity, 
pump photons with frequency ui p can downconvert into an entangled pair of photons, 
so long as the new photons' frequencies satisfy an energy conservation constraint, 

Ldp = UJ m + U„ , (1) 

where uj m and io n are the m th and n th resonant frequencies of the cavity. The range 
over which the downconvcrsion can occur is limited by the phasematching bandwidth 
of the cavity. 

More generally, multiple frequencies can be simultaneously phasematched inside 
the cavity, and upconversion is allowed as well as downconversion. The total 
Hamiltonian is then idealized by the following 

n = ihK^2 M -»»(«ln4-vJ. ( 2 ) 

p£P m+n=p 

where P is the (discrete) spectrum of the polychromatic pump, k is a global coupling 
strength, a\ is the creation operator for the n th mode, and M mn are matrix elements of 
a symmetric matrix. In units of the fundamental frequency of the cavity, the restriction 
of the inner sum to terms with m + n = p enforces photon energy conservation, and 
the magnitude of M mn determines the strength of the coupling in units of Hk, while 
the sign determines whether the photons are upconverted or downconvert ed. If all of 
the coupling strengths are equal, then M is a Hankel matrix whose elements are either 
or ±c, where c is a constant. Experimental simplicity further demands that there 
is no single-mode squeezing and hence we require that the elements along the main 
diagonal of M are all zero, i.e. M must be a hollow matrix. 

The Hamiltonian (2) can be used to create continuous-variable cluster states, 
defined as Gaussian states satisfying the relation 

P - Aq , (3) 

where p and q are vectors of amplitude and phase quadratures for the fc th cavity 
mode, qt- = at + a k and pk = — A is the symmetric weighted adjacency 
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matrix of the graph of the CVCS, and the arrow denotes the limit of large squeezing. 
Reference [8] first proved a nonunique relationship between the adjacency matrix A, 
which describes the cluster state, and the adjacency matrix M, which describes the 
coupling between the modes inside the cavity. The general relationship is somewhat 
complicated, so references [9, 7, 10] introduced the simplifying ansatz that 



In short, assuming A is an orthogonal matrix implies the graph describing the 
couplings between photons in the cavity (M) is identical to the graph describing 
the cluster state (A), which is a vast simplification over the general relationship. The 
derivation requires that M be the adjacency matrix of a bipartite graph, but as we 
will see, one can derive this from the Hankel and hollow constraints. 

We now see how the twin demands of experimental and theoretical simplicity 
for creating CVCS in an optical cavity naturally lead us to consider Hankel hollow 
orthogonal matrices, all of whose non-zero elements are either ±c, for some constant c. 
We can of course consider renormalized matrices instead by dividing out c, so that 
all entries arc ±1; after finding these renormalized matrices, we can then reintroduce 
the constant afterwards so that (4) holds. Thus, we see that finding and classifying 
Hankel hollow weighing matrices is our primary interest. 

For completeness, wc mention one further elaboration that is possible for this 
scheme. If the photons' polarization, transverse, or spatial degrees of freedom are 
taken into account, then the Hamiltonian (2) can be modified by adding additional 
indices for these modes. If these additional modes arc all frequency-degenerate, then 
an additional sum over these degenerate modes appears in (2), which allows for 
the symmetric intercoupling between the degenerate degrees of freedom. If such an 
interaction could be simultaneously phasematched by the nonlinear medium inside 
the cavity, then one is naturally lead to consider Wocfc-Hankel matrices, where the size 
of the blocks is equal to the number of degenerate degrees of freedom. We leave the 
consideration of block-Hankcl hollow weighing matrices to future work. 

3. Anticirculant graphs 

In this section we introduce the definition of anticirculant graphs and study some 
of their first properties. For the sake of self-containcdness, we will define here all 
the graph-theoretic concepts we require. For further background on theory of graphs 
the interested reader is referred to the book by Diestel [11]. Our reference about 
permutations and finite groups is Cameron [12]. 

We will work with simple graphs. A (simple) graph G = (V, E) is an ordered 
pair of sets defined as follows: V(G) is a non-empty set, whose elements are called 
vertices; E(G) is a non-empty set of unordered pairs of vertices, whose elements are 
called edges. An edge of the form {i, i} is called a loop. Two vertices i and j are said 
to be adjacent if {i, j} £ E{G). Often, this is simply denoted by writing ij. A graph 
is loopless if it is has no loops. The adjacency matrix of a graph G is denoted by A(G) 
and defined by 




(4) 



M = A. 



(5) 
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An adjacency matrix of a looplcss graph is said to be hollow. The degree of a vertex i in 
a graph G is the number of edges incident with the vertex, that is, #{j : ij £ E(G)}. 
A graph is regular if each of its vertices has the same degree. We say d-regular to 
specify that each vertex has degree d. A square matrix is Hankel if it has constant 
skew-diagonals. These are diagonals that traverse the matrix from North-East to 
South- West. (For this reason, these are also called antidiagonals). Notice that an 
n x n matrix has 2n — 1 skew-diagonals. An n x n Hankel matrix is said to be 
anticirculant (or, cquivalcntly, skew-circulant or backcirculant) if the i-th and (i + n)- 
th skew diagonals are equal. A permutation of length n is a bijection 7r : [n] — > [n], 
where [n] = {1,2, ...,n}. A permutation matrix of dimension n is a n x n matrix P 
with the following two properties: [P]ij £ {0,1} for every i, j; P has a unique 1 in 
each row and in each column. A permutation n is said to induce a permutation matrix 
P if [P]i l7r (i) = 1 for every i. 

Two graphs G and H are isomorphic if there is a permutation matrix P such that 
-4(G) = PA(H)P~ 1 . We write G = H to denote that graphs G and if are isomorphic. 
A graph G is bipartite if there is a permutation matrix P such that 

PA(G)P- = ( M ° T (6) 

for some matrix M. If A/ = J„, where J n is the all-ones matrix of size n, then 
G = K n>n , the complete bipartite graph. A graph is said to be anticirculant if it is 
isomorphic to a graph whose adjacency matrix is anticirculant. Let us denote by 
and by C n the set of all d-regular anticirculant graphs on n vertices and the set of 
anticirculant graphs on n vertices, respectively. Notice that #£ ra ,d = Q) ■ Therefore, 

En— 1 ^ — 1 / n\ 

Given G £ C n , the ordered set S(G) = (s : [4(G)] i. s = 1) is said to be the symbol of 
G, in analogy with the terminology used for circulant graphs. Indeed, S(G) specifics 
G completely. In cycle notation, each element s € S(G) is associated to a permutation 
of the form 

7T S = (l > s)(2 J s-l)...(a+l,n)(a + 2,n-l)--., (7) 

where the subtraction is modulo n. In line notation, tt s = (s)(s — l)(s — 2) • • • (n)(n — 
1) • • • (s + 1). On the basis of the definitions, we can collect the facts below: 

Proposition 1 Let G £ C n be a graph. The following statements hold true: 

(i) The graph G is d-regular with d = #S(G). 

(ii) If G is loopless then n is even. 

(Hi) If G is loopless then each s £ S(G) is even, 

(iv) The graph G is bipartite. 

Proof. In order: 

(i) Let 7r s be the permutation associated to the element s £ S(G). Then n s 
induces a symmetric permutation matrix P s . Since each P s is defined by s, 
we have [-P s ]-t,,j = 1 if and only if [Pt]i.j = for all t ^ s. Consequently, 
4(G) = J2seS(G) P" an d G is a regular graph with degree d = #5. 
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(ii) If the size of P s is odd, then [P s ]i,i = 1 for some i. This is because ir s is necessarily 
an involution of odd length, and for this reason it has at least one fixed point. 
(Recall that an involution is a permutation that is its own inverse.) 

(iii) Suppose on the contrary that s = 2k + 1, for some k > 1. Then, L4(G)] ls = 
[j4(G)]i 2fc+i = [A(G)] k+i,k+i = 1- contradicting the hypothesis of G being 
loopless. 

(iv) By the form of each tt s (see Eq. 7) , it is sufficient to observe that vertices labeled 
by odd numbers are adjacent only to vertices labeled by even numbers and viz. 



At this stage, we are ready to prove that the members of C n are Cayley graphs. 
The Cayley graph G = G(T, T) of a group T w.r.t. the set T C V is the graph in which 
V(G) = {T} and {g,h} G E(G), if there is s e T such that gs = h. The dihedral 
group D 2 k of order 2k is the group of symmetries of a regular /c-gon. The dihedral 
group D 2 k is nonabelian and presented as 

D 2k = (s, t : t k = e, s 2 = e, sts = t" 1 ), 

where t is a rotation and s is a reflection. The elements of D 2k are k rotations 
t a ,t,t 2 7 ...,t k ~ 1 and k reflections s,st,st 2 , ...,st . A graph G is said to be connected 
if there is no permutation matrix P such that PA(G)P^ 1 = ®^ Mj, for some matrices 
Mi. If a Cayley graph of a group T is not connected then it is the disjoint union 
of isomorphic Cayley graphs of a subgroup of T, with each isomorphic component 
corresponding to a coset of the subgroup. 

Proposition 2 Let G be a loopless anticirculant graph. Then G is a Cayley graph of 
the dihedral group with respect to a set of reflections. 

Proof. The adjacency matrix of a Cayley graph G = G(L, T) can be written as 
A(G) = J2t£T Preg(t); where p rcg is the regular permutation representation of T. 
When the order of T is n, this is an homomorphism of the form p TCg : T — > S„, 
where S n is the set of all n x n permutation matrices. Each p reg (t) describes the 
action of the group element t on the set {1, 2, n}. Let G <E C n be a loopless graph. 
By Proposition 1, we need to take n = 2k. Let us label the lines (rows and columns) 
of A(G) with the group elements of D n , of order n = 2k, in the following order: 

t° = e s t st t 2 st 2 ■■■ t k - x st k - x , 

where e is the identity element of the group. Consider p lcg (st l ), where 1 < I < k. By 
applying the generating relations of D n , it is straighforward to verify that p rog (s) = P 2 , 
and p reg (st l ) = P2+21, when 1 < I < k — 1. The elements of S(G) are then associated 
to reflections in D 2 k- 

2 — ► s 4 — ► st 6 — > st 2 ■■■ 2k — > st k ~ l . 

This concludes the proof of the proposition. ■ 

An example is useful to clarify this result. Let Z„ be the additive group of integers 
modulo n. The Cayley graph A(Z„, {1, n — 1}) is also called the n-cycle. This is the 
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unique connected 2-regular graph up to isomorphism. Let G G Cq^ be the graph with 
adjacency matrix 

/ 1 1 \ 

10 10 

1 1 

1 1 
1 1 

\ 1 1 j 

where s = (1, 2) (3, 6) (4, 5) and t = (1, 4) (2, 3)(5, 6). Then s 2 



A(G) = Prcg(s) + Prog (st) = 



= st 2 = e and 



(s, st) = Dq. Additionally, it is immediate to see that G = X(7,q, {1, 5}). 



4. Anticirculant weighing matrices 



In this section we will highlight the interplay between weighing matrices and 
anticirculant graphs supporting orthogonal matrices. 

A Weighing matrix W of order n and weight k, denoted by Win, k), is a square 
n x n matrix with entries [W]ij G { — 1,0,1}, satisfying WW T = kl, where I is the 
identity matrix. When k = n. then W(n, k) is said to be a Hadamard matrix. It 
is generally recognized that weighing matrices were first discussed by Frank Yates in 
1935 [13], while Hadamard matrices were introduced by James Sylvester and Jacques 
Hadamard during the second half of the nineteenth century (see [14]). Geramita and 
Seberry [15], and Koukouvinos and Seberry [16] are general surveys on this topic, as 
well as applications. 

We are specifically interested in circulant weighing matrices. A circulant weighing 
matrix of order n and weight k is a W(n, k) which is also a circulant matrix. In 
analogy with the previous section, we define the ordered set T(W(n, k)) = {s : 
[A(CW(n, k))]i, s = 1} to be the symbol of CW(n, k). Arasu and Seberry overview 
the subject in [17]. 

We shall give particular importance to graphs associated to weighing matrices. 
The graph of a (real symmetric) matrix M, denoted by G(M), is the graph defined 

by 



[A(G(M))]i 



if [M] itj 
if Mid 



7^0; 
= 0. 



The graph G(CW(n, k)) is a Cayley graph of the cyclic group Z ra with respect to the 
set T(W(n,k)). For instance, let us look at the matrix 

/ 1 1 1 -1 \ 



CW(6,4) 



V 




-1 
1 


1 
1 



1 



-1 
1 



1 



1 
1 



-1 
1 







1 
1 



-1 
1 



-1 
1 



1 
1 





(8) 



The graph G(CW(6, 4)) is the Cayley graph of Zq with respect to the set {1, 2, 4, 5}. 
This is illustrated in Figure 1. 

The following statement is one of the main tools in our discussion. For the 
purposes of this paper, it is sufficient to focus on matrices of even order. 
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Proposition 3 For each circulant weighing matrix M = CW(n, k) of order n = 21, 
with I > 1, there is a weighing matrix U such that G(U) £ H n ,k- The graph G(U) is 
loopless if each element of T(W{n,k)) is odd. 

Proof. Consider the permutation 

7r=(l,n) > (2,n-l),...,(M+l) (9) 

on a set [nj. Let be the permutation matrix induced by n. Labeling the rows and 
the columns of P n with the elements of [n] in the lexicographic order, we can write 

/ ■■• 1 \ 





V i 











1 



0/ 



Now let U = MP V , Since P^ and M are orthogonal matrices, U is also orthogonal. 
By the action of the permutation ir and the fact that M is circulant, the graph 
G(U) e £„, fc and S(G(U)) = {s = n - [A(CW(n, fc))] M + 1}, where t € T(W(n, k)). 
Notice that G(U) can have self-loops if we do not impose further restrictions. On 
the other hand G(U) is loopless if each s G S(G(U)) is even (Proposition 1). Such a 
property is satisfied only if each element of T(W(n, k)) is odd. ■ 

It is important to remark that P^ acts on columns and hence it does not generally 
preserve isomorphism of graphs. Proposition 3 prompts us to some definitions. An 
anticirculant weighing matrix of order n and weight k, denoted by AW(n, k), is a 
W(n, k) which is also anticirculant. A matrix AW(n, k) is a special kind of Hankcl 
weighing matrix. A Hankel weighing matrix of order n and weight k, denoted by 
M = HW(n,k), is a W(n,k) such that G(M) is anticirculant. While it is obvious 
that there exists an AW(n, k) if and only if there exists a CW(n, k), we need to take 
into account the following counterexample: 

Proposition 4 Not every HW(n,k) is an AW(n,k). 

Proof. Let P^ be the permutation matrix as in the proof of Proposition 3. The 
matrix 

P P 

P — P 



U 



(10) 
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is a HW(n, 2), for every n, since 

2Pl 



UU T = U 2 
because P^ = Pi . ■ 



IP' 



2 





' 21 


° ) 


H 


< o 


21 J 



Proposition 4 seems to indicate that classification theorems for circulant weighing 
matrices help only partially when attempting to classify Hankel weighing matrices. 

In the next two sections, our classification begins by considering graphs of small 
degree. This corresponds physically to the number of pump beams required to build a 
given CVCS. If a graph is r-regular, it requires 2r — 1 different pump beams to build 
the associated CVCS, since each of the 2r — 1 bands of the adjacency matrix requires 
a different pump frequency. Thus, we are primarily interested in graphs that have an 
interesting topology, but also have small values of r, since greater values would have 
greater experimental complexity. Cases of small r are also more tractable theoretically, 
which is another reason to focus on them. 



5. 2-regular graphs 

The unique - up to isomorphism - connected 2-regular graph on n vertices is the 
n-cycle, C n . A graph G is said to be the disjoint union of graphs Hi, ...,Hi if there 
is a permutation matrix P such that PA(G)P~ 1 = ® - =1 A(Hi). In such a case, we 
write G = 1+|' =1 H t . Equivalently, let G and H be graph such that V(G) D V(H) = 0. 
The set of vertices of G W H is V{G) U V(H) and the set of edges, E(G) U E(G). All 
useful information about graphs of matrices HW(n, 2) is stated as follows: 

Proposition 5 Let D n be the dihedral group of order n = 2k. Let T = {st l , sP} C 
D n , where 1 < i,j < k, be a set of reflections. There is a unitary matrix U such that 
X(D n ,T) = G(U) if and only if {i — j) mod k = (j — i) mod k. Moreover, if this is 
the case then k is even and therefore n is a multiple of A. It follows that the graph 
G G C n ^2 and G is the disjoint union ofk/2 copies of 'C4. 

Proof. It is known that given a group T and a set T = {g, h} C T, there is a 
unitary matrix U such that A(X(T,T)) = A(G(U)) if and only if gh^ 1 = g~ 1 h 
(see, e.g., [18, 19]). From this, sV-^sP)^ 1 = (st l )~ 1 sP. Since st l is a reflection, 
for every 1 < I < k, we can write (st 1 )^ 1 — sr. Thus sPsP = sPst 1 , that is, 
st l and sP commute. Since st l s = t k ~ l , for every 1 < I < k, it follows that 
t (k-i+j) mod k = t (k-j+i) mod k _ -phis proves that (i - j) mod k = (J - i) mod k, 
for st l and sP . Now, without loss of generality, assume that j > i. Note that 
(i — j) mod k = k — j + i. So, we need that j — i = k — j + i, which implies j — i = k/2. 
It follows that U has the form described in Eq. (10), up to isomorphism. ■ 

Figure 2 illustrates the graph X(D$, {st, st 3 }) associated to the matrix HW(8, 2). 
Since in general only a disjoint union of squares can be achieved with 2-regular graphs, 
and these graphs have only limited interest from a quantum computing perspective, 
we need to consider larger graphs. A proposal for creating such graphs was given in 
Ref. [9]. 
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Figure 2. The graph X(D S , {st, st 3 }) S C 4 tti C 4 . 



6. 4-regular graphs 

In this section we consider several constructions of 4-regular graphs of anticirculant 
weighing matrices. Unfortunately, we are unable to prove that this set of examples is 
exhaustive, so a complete characterization of the degree 4 case remains open. 

The CVCS graphs in this section can all be implemented with an experimental 
setup requiring no more than 7 pump beams. 



6.1. Graphs from AW {A, 4) 



4,4 



AW [A, A) = 



1 \ 



1 J 



The matrix 

/111 

1 1-1 
1-1 1 

V -l i i 

is a Hadamard matrix of order 4. Since the entries of a Hadamard matrix are all 
nonzero, G(M 4j 4) = K± , where denotes the complete graph on n vertices with 
a self-loop at each vertex. We can construct other anticirculant weighing matrices of 
higher order from M4 4. Given that we are interested in anticirculant matrices whose 
graph is without self-loops, the first row of these matrices must have nonzero entries 
only at even positions. Let ^4,4 be a family of matrices obtained from M^i be adding 
new lines (i.e., rows and columns) between the rows and the columns of M 4j 4. If we 
add the same number of lines between each column (resp. row), we obtain matrices 
of order Ak, for k > 1. For small orders, the first rows of these matrices are 

0,1,0,1,0,1,0,-1 

0,0,1,0,0,1,0,0,1,0,0,-1 

0,0,0,1,0,0,0,1,0,0,0,1,0,0,0,-1 

0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, -1 

Clearly, for each matrix M constructed in this way, we have G(M) G L n . When 
n = 4k, with k odd, the graph G(M) has self-loops. For this reason, we consider only 
the case n = Ak, with k even. In other words, we consider the case n = 81, for any I > 1. 
Let S(M) = {si, S2, S3, S4}. By the above construction s.; + i — Si = d+1 for i = 1,2, 3, 
and n — S4 + 2si = S2- This last condition provides that the number of columns 
between S4 and si is the same as one between S2 and s\ when looking at the matrix as 
wrapped on a torus, by gluing together the first and the last column. In this way, we 
can define a unitary matrix U such that G(U) = G(M). Let us denote by U\, U n 
the rows of U. The nonzero entries in U\ are exactly Ui fSl , U\_ S2 , Ui_ S3 and U\_ Si . In 
particular, C/i jSl = U\. S2 = Ui. S3 = 1 and Ui tSi = —1. Suppose s^+i — Si — 1 = d for 
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= S2- Thus the inner product between the rows U\ and 

U 1>Sl U hsi + U hS2 Ui, S2 + U hS3 Ui tS3 + U 1>ai Ui tS4 
U hSl U hS2 + Ui, S3 Ui, a3 + U Xj83 U hSi + Ui, SA Ui, ai 

1+1-1-1 

0. 

The distance d guarantees that each row (resp. column) of U has exactly 4 nonzero 
entries contributing to the zero inner product with other 3 rows (resp. columns). 
Orthogonality is guaranteed since no nonzero entries contribute to the inner product 
with all remaining rows (resp. columns). The numbers s\,...,s^ must be even and 
d = (n — 4)/4 must be odd, since n = 81 = (4d + 4). The matrices of order 8/ 
constructed in with this method will be denoted by M44 j. Essentially, orthogonality 
arises since we interlace M4 4 with itself d + 1 times. This can be seen directly in the 
matrix 

M 4A1 = AW(8,4) = 






1 





1 





1 





-1 \ 




























1 





1 





-1 





1 




























1 





-1 





1 





1 




























-1 





1 





1 





1 


-1 1 




















/ 



(11) 



which is constructed by interlacing two copies of M^. The boxed numbers are the 
nonzero entries in one of the two copies. At this stage, we can ask information about 
the graphs of matrices in M44,. It is sufficient to observe that 4 2 = M44 ® Fi + \, 
where the matrix 

/ •■■ 1 \ 



: 
V 1 ••• ) 

is 71 x n. Note that this is just the matrix from Proposition 3, but we have 
changed notation to make the dependence on the size of the matrix explicit. The 
matrix F n is the adjacency matrix of the disjoint union of n/2 graphs K<z, if n is even 
and the disjoin union of (n — l)/2 graphs K2 and a single vertex with a self-loop, 
if n is odd. The tensor product G ® H of graphs G and H is the graph such that 
A(G ® H) = A{G) ® A(H). The set of vertices of G ® H is the Cartesian product 
V(G) x V(H) and two vertices {u, u'} and {v, v'} are adjacent in G (8) H if and only if 
{u, v} & E(G) and {u 1 , v'} E E(H). The graph G ® H is connected if and only if both 
factors are connected and at least one factor is nonbipartitc. Here, K2 and K4 are 
both connected and K2 is bipartite. When I > 1, one of the factors is not connected. 
Therefore G^A^^j) is connected only if I = 1. In fact G(M4^ t i) = 1+J^ ; ^4,4, where 
K n n denotes the complete bipartite graph on 2n vertices. In Figure 3 is illustrated 
G(M 4 ,4,i); in Figure 4, G{M iAA ). 



all 1 = 1,2,3 and n — S4 + 2si 

Us2— SI ^ 

(U U U S2 - S1 ) = 
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Figure 3. The graph G(M 4j4 ,i) ^4,4 



Figure 4. The graph G(M 4i 4, 4 ) K 4ii ttl K 4 ^. 



The graph G{M i ^ i) is made of Z connected components. Can we find two 
permutation matrices P and Q such that G(PM^4jQ) ^ G(M^±i) and it is 
connected? In particular, for the moment, we consider only permutation matrices 
that preserves the distance d between the elements of S(M). From Ma^i, we can 
construct exactly I — 1 further matrices, by permutation the lines under this constraint. 
For example, when n = 24 = 8 • 3, the possible first rows are 

0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, -1 
0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, -1, 0, . 
0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, -1, 0, 0, 0, 

The position of the entry —1 does not matter, because we are interested in the graph 
G(-M4.4 j ;). Given s±, S4, in matrices of the form M^^j, when I > 1, 

gcd({s,; — Sj : for all i, j such that j > i} U 81) > 2. 

This implies that G(PM4 } 4jQ) = 1+J_^ ( A'4 4, for any two permutations P and Q 
preserving the distance d. 

6.2. Graphs from AW '(6, 4) 

It is natural to analyse the construction in the previous section, but replacing the 
matrix AVF(4,4) with some other anticirculant weighing matrix. Recalling that there 
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Figure 5. The graph G(M 6j4 ). 



Figure 6. The graph G(M 6:4 ,i) ^ G(M 6 , 4 ) (8) K 2 - 



is no j4W(5,4), the smallest available matrix of order n > 4 is AW(6,4): 



6,4 



>W(6,4) 



/ -1 
1 



1 

1 





V 



1 1 1 \ 
110-1 
110-11 
10-110 
0-1101 
-110 11/ 

illustrated in Figure 



This is obtained from the CW(6,4) in (8). The graph G(M 6j4 ) 
5, has 4 self-loops. 

The graph G(Mq^j), if loopless, is on \2l vertices, with I > 1. The graph 
G(-M6,4,i) = G(Me > 4) (£> K2 is connected since K2 and G(Mg,4) are connected and 
G(A/6,4) is nonbipartite. (The chromatic number of G(Mq^) is 3.) A drawing of 
G(M 6i 4,i) is in Figure 6. 

In general, Me,4,i = Mqa <g> F1+1. If we avoid self-loops, 



G(M M ,0 = |±| j G(M M ,i)- 
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By permuting the lines of M g ,4 (shifting to the right the columns), we have the matrix 



/ 



ML 



AW(4,4)i = 



V 



-1 1 



/ 



Apply our construction to M\ 4 , we obtain graphs on 6k vertices, with k odd. However, 
since G(M 6i4j i) = G(Afg 41 ), these graphs will be a disjoint union of copies of 
Gr(M6,4,i). So the permutation is of no use if our intention is to get different graphs 
from M6,4- On the other side, the structure of M\ A suggests a different construction, 
which will allow us to obtain matrices of the form AW (Am, 4), whose graph is loopless 
and connected for every m > 3. Notice that 



where 



M M = 



.4 



.4 
D 



D 
A 




B = 




Then 2AB = and A 2 + B 2 



M 6A M 6A 



41 and indeed 
A 2 



(M 6A ) 



^B 2 
2AB 



2AB 

A 2 + B 2 



41. 



We can add an even number of extra lines to Mq^ and still preserve its block structure. 
For example, the matrices 



A' 



( 





-1 





1 





^ 




( 


1 





1 










-1 





1 













1 





1 
















1 











1 







1 











-1 




1 











1 





, Bf = 


1 











-1 
















1 





1 













-1 





1 


\ 








1 





1 


J 




I 





-1 





1 






satisfy the above conditions. Hence the matrix 



Mi 



A' B' 
B' A' 



is an AW(12,4). In general, let Ai6,4 be the family of matrices obtained with this 
construction. If M G Mq A then s\ = 2, s-i = 4, S3 = n/2 + 2 and S4 = n/2 + 4, 
where S(M) = {si, s 2 , s 3 , S4}. These are AW(4k, 4) when s 3 — s 2 = (n — 6) /2 is an 
odd number. The construction gives also matrices of order 21 > 6, for any I. Yet the 
associated graphs will have loops when I is odd. This is the reason for taking multiples 
of 4. If M G A46.4 then G(M) is connected, in virtue of the fact that s 2 — si = 4— 2 = 2 
and then gcd({s,; — Sj : for all i,j such that j > i] U 4k} = 2. 

The graphs of matrices in the set Mq A , with order 4m, are the direct product of 
Cim with the graph . (Recall that is the complete graph on n vertices with 
self-loops; these are the graphs with adjacency matrix J„.) The direct product G x H 
of graphs G and H has set of vertices V(G x H) = V(G) x V(H) and two vertices 
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Figure 7. The graph G(M 7j4 ). 



Figure 8. The graph G(M 7 ,4,l)- 



{u, u'} and {v, v'} are adjacent in G x H when {it, v} € E(G) and {V, v'} £ E(H). By 
the definition, A(G x iJ) = A(G) 8) For this reason G x iJ is sometimes called 

the Kronecker product of graphs (see Imrich and Klavzar [20], Ch. 5). If G G -Mqa 
then G = K} x G 2m and M(G) = J 2 ® A(C 2m ). 



6.3. Graphs from AW (7, 4) 
Let us consider the matrix 



M 7j4 = AW{7,4) = 



\ 



1 

-1 


1 


1 



-1 




1 



1 
1 



The graph of M-j^ is in Figure 7. 

Interlacing Mj t 4 with itself, by the same method used for (see section 6.1, 

Eq. (11)), we obtain the matrix 

M 7AA = AW(U,4) = AW (7, 4) ® F 2 (12) 
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Figure 9. Hamiltonian cycles in graphs of matrices AW (12, A) and 7iW(14,4). 



The graph of Mf^i is in Figure 8. 

This is on 14 vertices, 4-regular and connected. At a first analysis, G{M-r^\) 
does not seem to be a direct product of graphs. It is nontrivial to extend G(My^ j i) 
to an infinite family of connected regular graphs, as we did by taking Mq^. Examples 
show that graphs obtained by Mj^ are not connected, but G(M- 7y 4^). 

7. Open problems 

• Study the structure of the graph G{M-j^{). Determine if it belongs to a family of 
connected graphs arising from AW(7, 4), having an infinite number of members, 
as it is for the graphs arising from AVT(6, 4). 

• Since there is no general classification of weighing matrices, characterizing graphs 
of Hankel weighing matrices seems to be out of reach. This could be however 
done for graphs of degree 4 and 9, on the basis of present knowledge on CW(n, 4) 
and CW{n,9) (see [17]). 

• We have seen that there are examples of Hankel weighing matrices which are 
not anticirculant. A way to study the relation between these matrices and 
anticirculant ones, would be to prove that the graphs of Hankel weighing matrices 
are always isomorphic to graphs of AW(n, d), except for some special cases which 
will include the examples in proposition 4. 

• A Hamiltonian cycle in a graph is an ordered set of sequentially adjacent vertices, 
in which every vertex of the graph appears exactly once. Graphs with Hamiltonian 
cycles are said to be Hamiltonian. Alspach and Zhang [21] conjectured that 
a connected Cayley graph G(D n ,T) is Hamiltonian whenever T is a set of 
reflections. According to this conjecture, any connected anticirculant graph is 
Hamiltonian whenever it is looplcss. Additionally, Gutin et al. [22] conjectured 
that graphs of unitary matrices arc Hamiltonian if connected. According to these 
conjectures the graphs of Hankel weighing matrices arc Hamiltonian if connected. 
This is the case for the graphs considered here. In Figure 9 are drawn the graphs 
K% x C$ and G{Mj^^). The highlighted edges represent Hamiltonian cycles. 
Perhaps the physical picture that motivates our study could provide some insight 
toward proving these conjectures. 

• We have seen that graphs of the form x Ci m are graphs of unitary matrices. 
The quantum dynamics induced by these matrices is trivial because of symmetry. 
However symmetry can be broken with the action of a diagonal matrix without 
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altering the zero pattern of the unitary. This suggests the possibility of 
constructing unitaries with a nontrivial dynamics whose graph is a undirected 
Cayley graph of the dihedral group. 

• Sections 5 and 6 make extensive use of the following trick. Begin with a circulant 
weighing matrix, then reverse the columns. If the resulting graph has loops, one 
can create a loopless graph by the "interlacing" trick. In fact the interlaced graph 
will again be connected if the original graph was connected but not bipartite. 
Can one classify the Hankel and anticirculant weighing matrices which cannot be 
obtained by this trick? 

8. Conclusion 

We have demonstrated connections between an approach to quantum computation 
with optical modes and the existence of certain weighing matrices with Hankel 
structure. Because the degree of the graph associated to a matrix corresponds to the 
experimental resources required to implement the graph as a CVCS, it is important 
to characterize which graphs with small degree have adjacency matrices with Hankel 
structure. We proved some general theorems about such matrices, classified completely 
the case of degree 2, and provided some examples with degree 4. We also raised some 
interesting open problems that might help shed some light on which graphs can be 
implemented in a single OPO using the scheme in Ref. [7, 10]. 
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